Localized vibrational modes in optically bound structures 
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We show, through analytical theory and rigorous numerical calculations, that optical binding can organize 
a collection of particles into stable one-dimensional lattice. This lattice, as well as other optically-bound 
structures, are shown to exhibit spatially localized vibrational eigenmodes. The origin of localization here 
is distinct from the usual mechanisms such as disorder, defect, or nonlinearity, but is a consequence of the 
long-ranged nature of optical binding. For an array of particles trapped by an interference pattern, the stable 
configuration is often dictated by the external light source, but our calculation revealed that inter-particle 
optical binding forces can have a profound influence on the dynamics. © 2008 Optical Society of America 
OCIS codes: 140.7010,220.4610,220.4880,999. 9999(Optical Binding) 
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Since its introduction many years ago,^ optical manip- 
ulation has evolved into a major technique for manipu- 
lating small particles, and recently, simultaneous manip- 
ulations of multi-particles have been demonstrated. ^It 
is known that in addition to the well-known one-body 
force such as the gradient force that depends on the in- 
tensity profile, there is an optical binding (OB) force 
that couples the particles together. •^■'^ Nevertheless, for 
an extended array of particles, the nature of OB is not 
fully understood, although some theoretical efforts were 
devoted to small clusters.^ '^ As the principles underlin- 
ing these intcr-particlc forces are different from that of 
the traditional light-trapping, we expect some new and 
interesting applications. 

In this paper, we demonstrate an interesting conse- 
quence of OB in a spatially extended structure bound 
by light: the existence of spatially localized VEM (vibra- 
tional eigenmodes). We illustrate the physics by consid- 
ering a one-dimensional "lattice" bound by light. Wave 
localization is known to occur in defect or impurity sites 
of an otherwise ordered lattice. In solids, the "defect" 
can be impurity atoms that localize phonons, and in the 
intrinsic localized modes, the "defect" is derived from 
the nonlinearly excited particles.^ Here the localization 
occurs in the linear dynamics regime in an ordered array 
of identical particles without defect or disorder. 

Optically bound structures have been investigated in a 
number of recent experiments. Stable cluster configura- 
tions had been realized^'^'^'^'^^'^^ and vibrational mo- 
tions were observed.^ In particular, the most commonly 
observed geometry is an one-dimensional array of parti- 
cles, bound by a pair counterpropagating beams^'*'^ or 
evanescent waves. ^"^^^ 

Consider a linear chain of N evenly spaced spheres 
in air. The particles have mass density p=l,050 kg/m'^, 
dielectric constant £=2.53 (^^polystyrene) , and radii a = 
A/10 = 52 nm, so that they are small compare to the 
incident light's wavelength A=520 nm. The particles are 
illuminated by the standing wave formed by a pair of 



counterpropagating plane waves 



Er 



,{r) ~ 2Eo COS {kz)x, (1) 

where k is the wavenumber, and the intensity for each 
beam is set to be 0.01 W/^m^.'^''^ 

To calculate the optical force acting on the particles, 
we employ the rigorous and highly accurate multiple 
scattering and Maxwell stress tensor (MS-MST) formal- 
ism,^ which requires no approximation and subject only 
to numerical truncation errors (we use multipoles up to 
L=6). The optical force tends to drive small particles 
to the region of strong light intensity. For an array of 
N evenly spaced particles aligned along the z-axis, one 
expects a stable one-dimensional lattice with a lattice 
constant of A/2: 



Rr, 



(0,0,nA/2), 71 = 1,2..., TV 



(2) 



where i?„ is the equilibrium position for the n-th parti- 
cle. Indeed, we found that the geometry defined in Q 
corresponds to a zero-force configuration and the con- 
figuration is proven to be stable by using linear stabil- 
ity analysis.**' ^^ The longitudinal trapping (along the z- 
axis) is mainly provided by the gradient force of the inci- 
dent beam, and it is further enhanced by OB.^^ On the 
other hand, the transverse stability (on the x^z-plane) 
is solely induced by OB. We note that there are other 
beam configurations, other than that specified by ^, 
that can stabilize a linear chain as demonstrated by re- 
cent experiments. ^'®'^*'^'* 

The VEMs are obtained by diagonalizing the force 

matrix {K)jk = d{f i^ i^^)j/d{Ax)k,'^ which is found 
by linearizing the optical force near the equilibrium: 

flight ~ KAx, where Ax is the displacement vector 
of the i-th particle away from its equilibrium configura- 
tion. The vibration profile of the VEM is described by the 

eigenvectors V oi K, and the natural vibrational fre- 
quency is flQi = {—Ki/mY''^ where Ki is an eigenvalue 

of K and m is the mass of a sphere. Due to the reflection 
symmetry, the modes fall into three separate branches 



(each of N modes) , corresponding to the vibrations along 
the three Cartesian directions (see Fig. ^e)). We shah 

denote the branches as the fc-branch, iJ-branch, and 

(fc X i?)-branch, corresponding respectively to particle 

displacements along the incident wavevectorfc = ±fcl, 
the incident polarization (x-axis), and the y-axis. 

The degree of localization of the modes can be quan- 
tified by calculating the inverse participation ratio^^ 



To leading orders, the force matrices for the three 
branches, evaluated using the P.E. model, are 



(^^-.™,..)'^ 
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(3i\l^q\n)-' 






(5) 



(/.P.i?.), = 
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N 



( AX^*) Ar^*) AZ^*) 



(3) 

which indicates the number of particles participating the 
vibration. Here, the index i stands for the i-th eigenmode 
and AX,V is the vibration amplitude of the n-th parti- 
cle along the cc-axis. A small value of I.P.R. indicates 
a localized mode, while I.P.R. ~ N indicates a delo- 
calized mode. Fig. Q] shows the I.P.R. computed by the 
MS-MST formalism. For comparison, the I.P.R. for an 
ordinary "ball and spring" model is also plotted in Fig. 
n^d), where a lattice of 100 particles are connected to its 
nearest neighbors by a Hooke spring. As expected, the 
ball and spring model supports only propagating modes 
in which the displacement of the n-th particle ^ g™?^^ 
where q is the phonon wavevector and A is the lattice 
constant. Depending on whether gA is an integer multi- 
ple of TT, I.P.R. takes cither -200/3 or -100. 

In general, the VEMs of the optically-bound lattice 
are more localized than the propagating modes, espe- 
cially for the fc-branch. A few modes selected from the k- 
branch is shown in Fig. [5] The high-frequency modes are 
highly localized near the center of the lattice (e.g. Fig. 
El^c)), while those with a lower vibrational frequency are 
less localized (e.g. Fig. |2Jd)-(e)). For very low frequen- 
cies, the modes are further delocalized spatially (e.g. Fig. 
I^ff)), with the vibration being stronger on both ends. 
The evolution of a VEM as the number of particles in- 
creases is also depicted in Fig.|2Ia)-(c); clearly the overall 
profile of the modes are getting more and more localized 
as the number of particle increases. 

The physics of the localized mode (LM) can be cap- 
tured qualitatively by a simple potential energy model 
(P.E. model)/* For small (a ^ A) lossless dielectric par- 
ticles placed in a standing wave of light, one may define 
an approximate potential energy for the light-induced 
mechanical interaction as 



[/ = -£ (a/4)|i,;„(f„)p-aV2 



n— 1 m<Cn 



where a = 4TTeoa^{e — l)/(£ + 2) and 



G{R) = e''="/47reoi?' 
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and 




4(|Z-n|7r)"^ 

-9(|/-n|7r)"^ 

-9(|/-g|7r)-5] 



(1 = 1) 
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where 

(3 = A;^a^i?o/27reo, and I and q are particle indices. The 
I.P.R. computed using the P.E. model is plotted in Fig. 
Qas dotted lines, which are surely not quantitative com- 
pared with the exact result, but nevertheless captures 
the salient features of the rigorous calculations. 

It is evident from ^ and that the modes of the 

{k X _B)-branch and the _E-branch are similar, because 
the leading terms arc essentially an action-reaction cou- 
plings between every pair of particles, with the coupling 
strength being proportional to inverse-cubic distance. 
These two branches are more localized than those of 
the ball and spring model because the interaction has a 

longer range. *^ The fc-branch is the most localized and 
interesting. Its force matrix consists of two components, 
the long range (inverse distance) action-reaction cou- 
pling and Kiocai{l) which acts like a spring that ties the 
l-th particle to its equilibrium position. The first term of 
Kiocai{l) is caused by the incident beam and is the same 
for each particle. This term gives a frequency gap at low 
frequency (e.g. between and 4.7 MHz in Fig. ^a)), 
while the second term is induced by OB. One may de- 
fine an intrinsic vibration frequency for every individual 
particle as 



^ntrinsic(0 - \^ -Kiocai{l)/m, (8) 

plotted in Fig. Efg). We note that the first term of 
Kiocai{l) contributes a constant to f^jntrinsic^^)' ''^tiil*^ 
the term due to OB gives a position dependent con- 
tribution that makes f^jntrinsicC) higher (lower) near 
the center (ends) of the lattice. It is the variation of 
^intrinsic(0 along the chain that elicits the enhanced 



localization: only particles near the center (both ends) 
participate in the high (low) frequency vibrations, see 
Fig. lie) (Fig.iaf)). 

We now consider the strength of the OB. As revealed 
by recent theoretical^' ^^ and experimcntaP'^'^'^^'^^ 
works, the optical force on microspheres can dominate 
over other relevant interactions such as gravity, van der 
Waals, and thermal fluctuations. For the lattice consists 
of smaller spheres defined in ^, the potential energy 
per particle U/N for A^^l, 10, 50, 100 are respectively 
-9.6, -10.5, -11.3, -11.7 kBTnoom, and the chain should 
thus be thermally stable at the assumed intensity. Fur- 
thermore, U/N is enhanced by more than 20% as N is 
increased from 1 to 100, implying that the OB carries a 
non-negligible contribution. 

We have showed that OB can bind a collection of 
particles into a ID lattice that is stable in all three 
dimensions. We shall emphasize that the localization 
discussed here is a general phenomena for optically- 
bound structures that are spatially extended, and it 
is not restricted to the particular geometry or incident 
wave considered here. We found that LMs are also ob- 
served in other structures such as the photonic cluster 
made from microspheres shown in Fig. 4(g) of refer- 
ence 4, and also another lattice configuration defined 
by ^„ = (0,nA,0), n== 1,2..., iV ■ A difference be- 
tween this lattice configuration and (O is that the lattice 
constant of the later (former) is dictated by optical bind- 
ing (trapping). It is the long-ranged OB that induces the 
variation of f^intrinsic (0 ' "^^lich in turn induce the local- 
ization. 

It is worth to note that in the case of the ID array 
specified by (jSJ, the stable configuration is defined by 
optical traps produced by the incident wave rather than 
OB, yet OB plays a crucial role on the dynamics. The 
quasi-stable dynamics that arises from the nonconserva- 
tive nature of the optical forces,^ and the LMs consid- 
ered here, could be major causes of the inconsistencies 
between the vast amount of light-trapping experiments 
and theoretical predictions where OB is neglected.^ A 
deeper investigation into the subject would be an inter- 
esting and important research topic for the future. 
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Fig. 1. (Color online) Natural vibration frequencies fioi 
versus the inverse participation ratio, for the ID lattice 
with A^ = 100. Panel (a), (b), and (c) correspond re- 
spectively to fc-branch, [k x ii^)-branch, and i?-branch. 
The open circles are obtained by the rigorous MS-MST 
formalism and the dotted line is that of the P.E. model 
(gj. Panel (d) shows the I.P.R. for the "ball and spring" 
model. Panel (e) shows schematically the direction of the 
particles' displacements for the three branches. 
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Fig. 2. (Color online) The profiles of a few selected VEMs 

in the fc-branch. Dashed lines show the equilibrium po- 
sitions. Each circle represents one particle. Panel (a)-(c): 
the highest frequency mode for a lattice containing (a): 
Af = 10 particles, (b): A^ = 50, and (c) A^ = 100. Panel 
(c)-(f): the VEMs for A^ = 100, with (c) showing the 
highest frequency mode, and (d)-(c) correspond to two 
intermediate frequencies, and (f) shows the lowest fre- 
quency mode. Note that the intcrparticlc separation and 
the size of the particles are not drawn to scale. Panel (g) 
shows ^intrinsic '^ith A^ = 100, see (jSJ. 
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